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INTRODUCTION 


A computer program based on state-of-the-art compressor and 
structural technologies applied to bladed shrouded discs has been 
developed and made operational in NASTRAN Level 16. 

The problems encompassed include aeroelastic analyses, modes 
and flutter. 

The program is documented in the form of five NASA Contractor's 

Reports one Technical Report and four Updates to NASTRAN Level 

16 Theoretical, User's, Programmer's and Demonstration manuals. 

This report describes the Demonstration manual updates. 


DEMONSTRATION MANUAL UPDATES 


|P RIGID FORMAT No. 16* Static Acrothermoelast ic 1 PcsIgn/AnaVysIs 1 

l ♦ 

[ Aerot’astic 'Design/Analysis' of an Axial Flow Compressor Stage (16-1) 

f 

6 . ’ A* Description 

i 

~ _ This problem Illustrates the use of the static aeroelastlc analyses of the 

first stage rotor of an axial flow air compressor to determine, 

( 

I 1) the "as manufactured" blade shape required to produce the design point 

j pressure ratio ("design" problem), and 

! 11) the operating point of the "flexible" designed blade ("analysis" 

I 

problem). The total stiffness matrix, consisting of the elastic and geometric 
stiffness matrices, at any off-design operating point is saved for use In subse- 
k quent modal and flutter analyses. (See Demonstration Manual example 9-1), i 

( The 43-blade rotor is designed to develop a total pressure ratio of 1,85 at 

a speed of 16043 rpm and an air flow rate of 73.15 Ibm/sec,. The finite element 
model of a representative sector of the rotor Is shown In Figure 1, 

B. Input 

Bulk o.ata cards used include DTI, STRCAML1, PARAMeters APRESS, ATEMP, FXC00R, 
FYCflJlR. FK00R, 1PRTCF, IPRTCI, IPRTCt, KT0UT , PGE0M, SIGN, STREAML and 20R1GN aS 
Illustrated In the User's Manual Section 1*15.3. 

W 

C . Analyses and Results 

, The rigid blade of Figure 1 produces a total pressure ratio of 1.85 at 

r 

16043 rpm and 73.15 lbm/sec air flow rate (Table 1), Because of the elasticity 
of the material, and under the action of centrifugal and aerodynamic pressure *nd 
thermal loads, t:ie blade deforms and produces a total pressure r ) t i o greater than 
the design value. A "redesign" of the rigid blade, considering the elastic and 
geometric properties of the bladed disc sector, enables determination of the "as 
manufactured" blade shape that, when loaded and deformed, would produce the design 
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pressure ratio. The 'rigid' performance of the "as manufactured" blade shape 
obtained at the end of the Design problem is also shown in Table 1, 

This blade shape is then "Analyzed", in the current demonstration example, at 
the same (design) speed and flow rate to determine the 'flexible' operating pressure 
ratio. This value (1,84) can be further improved to approach the desired (1.85) 
pressure ratio by reducing the Parameters FXC00R, FYC00R and FZC00R in the Design 
problem (see User's Manual Section 1,15,3). 

The blade shape at various stages during the Design and Analysis problems, as 
reflected by the grid point coordinates, is also shown in Tables 1 and 2 . The 
coordinates are expressed in the basic system of figure 1, 
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Axis of Rotation 


Figure 1 



Finite Element Model of an Axial Flow Compressor 
Rotor Sector, and the Basic Coordinate System 


16 . 1-3 ( 9 / 30 / 78 ) 


Table 1 » Design Problem 


Total Pressure Ratio 
Rotational Speed, rpm 
Air Flow Rate, lbm/sec 


Grid Points 


1 

2 

3 

4 

5 

6 

7 

8 
9 


Ini tial , Designed Blade 
'Rigid' Performance 

"As Manufactured" Blade 
'Rigid' Performance 


1.85 



U80 



16043 

\ 


16043 



73,1 5 



73,15 


X, in. 

Y,1n. 

2 , in. 

X, In. J 

' Y, in. 

maam 

-0,8900 

-0.2732 

3,7902 

-0.8981 ; 

-0.2755 

3.7796 

0,0 

O.Q532 

3.9996 

-0.0001 

0.0540 

3.9986 

0.8980 

-0,2499 

4.1926 

0.8979 

-0.2464 


-C ’630 

-0.5004 

5.4772 

-0.7653 

-0.4830 

5,4554 

0,0 

0.0209 

. 5.5000 

-0.0005 

0.0209 

5,4985 

0.7800 

0,2342 

5,4950 

0.7799 

0.2307 

5.4889 

-0.6290 

-0,7494 

7.3620 

-0.6386 

-0.7217 

7.3281 

0.0 

0,0131 

7,4000 

-0.0091 

0.0155 

7.3976 

0.6290 

0.6369 

7.3725 

0,6240 

0.6123 

7.3416 

-0.4240 

-0.9330 

9,9564 

-0,4058 

-1.1351 

9.8905 

0.0 

-0.0221 

10.0000 

-0.0106 

-0,0236 

9.9970 

0.4240 

0.7598 

1 

9.9711 

0.4140 

0.8134 

9.9304 
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total Pressure Ratio 
Rotational Speed, rpm 
Air Plow Rate, 1 bn/sec 


Table 2, Analysts Problem 


"As Manufactured" Blade 
'Rtgt<:' Performance 

"As Manufactured" Blade 
'Flexible* Performance 


1.80 



1.84 



16043 



16043 



73,15 



73,15 


X » 1 n , 

Y,1n, 

Z,1n, 

X , i ri , 

Y , t n , 

Z,tn. 

-0.8981 

-0,2755 

3.7796 

-0.0979 

-0,2814 

3.7712 

-0.0001 

0.0540 

3.9986 

0.0001 

0.0516 

4.0003 

0.8979 

-Q.2464 

4.1847 

0.8901 

-0.2461 

4.1795 

-0.7653 

-0.4830 

5.4554 

-0.7726 

-0.4.744 

5,4413 

-0,0005 

0,0209 

5.4985 

-0,0031 

0.0228 

5.5033 

0.7799 

0.2307 

5.4889 

0.7797 

0.2247 

5.4889 

-0,6386 

-0.7217 

7.328) 

-0.6646 

-0.7082 

7.3062 

T 

O 

o 

o 

vD 

0. 0155 

7.3976 

-0,0157 

0.0164 

7,4058 

0.6240 

t, 6123 

7.3416 

0.6303 

0.5962 

7,3237 

-0.4058 

-1.1351 

9.890/{ 

-0,5237 

-1.1552 

9.8520 

-0.0106 

-0.0236 

9.9970 

-0.0320 

-0.0656 

10.0079 

0.4140 

0.8134 

9.9304 

0.4130 

0.7329 

9,9093 
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RIGID F0RMAT Ho, 9 (APP AERp), Aeroelastic Analysis 


Modal i Flutter and Subcrltical Roots Analyses of an AxUl Flow 

| 

Compressor Stage (9-5-1) 

A. Description 

The problem Illustrates the use of the aeroelastic cyclic modal and flutter 
analyses of the first stage rotor of an axial flow air compressor to* 

I) determine the natural frequencies and mode shapes of the bladed disc 
sector, of Figure 1, which exhibits rotational cyclic symmetry. The total stiff- 
ness matrix, including the differential stiffness effects at the operating point 
under consideration, saved during the Static Aerothermoel astic "Analysis" (see 
Demonstration Manual example 16-1) Is used for the cyclic modal analysis, 

II) examine If the operating point being considered Is a flutter point by 
analyzing the V-g and V-f plots basea on user-selected combinations of densities, 
inter-blade phase angles and reduced frequencies, and in the process 

III) Identify the subcrltical (stable) roots. 


B. Input 

Bulk data cards used include AERP, FLFACT , FLUTTER, MKAE&pl , S7REAML1 , 
STREAML2 and PARAMeters IREf, KGGIN, LMpDES, MAXMACH, M HIM AC II, MTYPE and PRINT as 
described in the User's Manual Sections 1.15.2 and 1.15.5, Bulk data cards 
CYJpIN and PARAMeters CTYPE, KINDEX and NSEGS are discussed in Section 1.12 
of the User's Manual. 

* 

C. Analyses and Results 

The finite element model of the bladed disc sector analyzed is shown in 
Figure 1. The first five zeroth harmonic natural frequencies and mode shapes of 
the sector are noted in Table 1. The grid points on the hub in contact with the 
compressor shaft were permitted radial translational degree of freedom only, 
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As » typical example, the first of the three frames of V-g and V-f plots out- 
put requested In this demonstration problem Is shown In Figure 2, The density 
and Interblade phase angle are held constant at (Q. 059 x 1,507 E-6) si Inch/in 3 
and 180°» respectively, for this frame. The three reduced frequencies are Identi- 
fied by the symbols *(k»0.3), C(k«0.7) and A(k*1.0), The flutter summary for the 
three (p,o, k) groups Is presented in Table 2, 

A close examination of the damping curves shows that the damping Is nearly 
zero In the fourth mode of frequency 1797 ll 2 . 

The Implied density, velocity and reduced frequency are, respectively 
(0,059 x 1,507 E-6) $linch/in 3 , 1.055 E4 in/sec and 1,0 as compared with the 
actual values of these quantities as (0.059 x 1,507 E-6) $11nch/in 3 , 1.910 E4 
In/sec and 1,0, respectively. 

The ratio t ed^^ac tua 1 not be * n 9 equal (or close) to 1,0 discounts the 

current operating po*nt as being on a flutter boundary, at which all the ’ihree 
Implied quantities equal the actual quantities. 

The demonstration example discussed has been presented principally for the 
purpose of illustrating the procedure for axial flow compressor flutter analysis. 
In order to locate the unstallcd flutter boundaries over the entire region of 
operation of the compressor stage, similar analysis would be required for a 
series of operating points, harmonic numbers, interblade phase angles and 
reduced frequencies for both the stage rotor and the stator. Appropriate super- 
position of the rotor and stator results, would then help Identify the unstallcd 
flutter boundaries on the compressor stage map. 
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Axis of Rotation 


Figure 1. Finite Element Model of an Axial Flow 
Compressor Rotor Sector* and the Basic 
Coordinate System 
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Figure 2. V-g, V-f Plots? Frame 





















